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ABSTRACT 

Gamma ray emission from dark matter subhalos in the Milky Way has long been sought as a sign 
of dark matter particle annihilation. So far, searches for gamma-ray continuum from subhalos have 
been unsuccessful, and line searches are difficult without prior knowledge of the line energies. Guided 
by recent claims of line emission at 111 and 129 GeV in the Galactic center, we examine the coadded 
gamma-ray spectrum of unassociated point sources in the Second Fermi-hAT catalog (2FGL) using 
3.9 years of LAT data. Using the SOURCE event class, we find evidence for lines at 111 GeV and 129 
GeV with a local significance of 3.3tr based on a conservative estimate of the background at E > 135 
GeV. Other 2FGL sources analyzed in the same way do not show line emission at 111 and 129 GeV. 
The line-emitting sources are mostly within 30 degrees of the Galactic plane, although this anisotropy 
may be a selection effect. If the double-line emission from these objects is confirmed with future data, 
it will provide compelling support for the hypothesis that the Galactic center line signal is indeed 
from dark matter annihilation. 

Subject headings: gamma rays — line emission — milky way — dark matter 



1. INTRODUCTION 

Although a wide variety of cosmological and astro- 
physical observations provide compelling evidence that 
nonbaryonic dark matter constitutes ~ 80% of the to- 
tal matter i n the Universe, w e still know little about its 
nature (e.g. Bergstrom 20121 : Hooper fc Profuniol 2007 ; 



Bertone et al.ll2005l ). 

In many models of weakly interacting massive parti- 
cle (WIMP) dark matter, particles annihilate and/or de- 
cay to gamma rays directly or indirectly. Gamma-ray 
photons at energies E 3> 1 GeV travel in straight lines 
without significant energy losses in the local Universe, 
allowing their spatial distribution to serve as a tracer 
of the dark matter distribution. Regions of high dark 
matter density such as the Galactic center, galaxy clus- 
ters, and dwarf galaxies have been suggested as possible 
sources. In addition, many DM subhalos in the MW 
may shine in gamma rays and have no counterpart at 
other wavelengths, making them promising sources (see 



involving one photon plus a Higgs boson, Z boson, or 
other chargeless non-SM particle. No plausible astro- 
physical background can pro duce a line, although a nar- 
row feature is possible (see lAharonian et alj|2012l ). In 
most models, the line flux is suppressed by a loop fac- 
tor relative to the contin uum, implying it is i 2-3 o rders 
of magnitude fainter (e.g. Bergstrom fc Ulliol 1997). Al- 



though this is not true in all models (e.g. [B ergstro m et al 
199a Bergstroml2000HBertone et al.ll2009Ujackson et al 



recently e.g. Bclikov et al. 


2011 


Ackermann et al.ll2012b 


Mirabal et al. 20121 and r 


eferences therein) 



The primary challenge in such searches is to under- 
stand the background from conventional astrophysics 
well enough to distinguish a dark matter signal. A 
"smoking-gun" signal of annihilating dark matter would 
be the discovery of one or more gamma-ray lines. The 
line(s) could be produced by dark matter decays or an- 
nihilations into two photons, or two-body final states 
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2010tlClinell2012l ). this theoretical prejudice led previous 
studies to focus on continuum searches. However, tenta- 
tive evidence for gamma-ray line emission at ^130 GeV 
toward the inner Galaxy has been fou nd with 3-3g s ig- 
nificance after trials factor 4 correction 

In our recent paper (|Su fc Finkbeinerll2012bl ). we have 
performed a study with various data analysis methods 
and obtained 6.5c local significance of the gamma-ray 
study with various data analysis methods and obtained 
6.6a local significance of the gamma-ray line structure, 
and 5.0-5.5(7 after trials factor (depending on whether we 
assume one line or two lines). In fact, we found two lines 
centered at 111 GeV and 129 GeV provide a better fit to 
the data. 

The high significance of this result does not address 
concerns about instrumental artifacts, such as energy 
mapping errors that could give rise to spectral bumps 
and dips (jFinkbeiner et al.1120121 ). Such concerns must 
be addressed by an independent analysis of photons from 
other parts of the sky. For example, detection of lines at 
111 and 129 GeV elsewhere on the sky in multiple unas- 
sociated LAT sources, but not in any class of associated 
LAT sources would rule out an energy mapping error in 

4 Also known as the "look elsewhere effect" 




Fig. 1. — Upper left: the spectrum of 16 unassociated point sources in the 2FGL catalog, chosen to have a photon at E = 100 — 140 
GeV (vertical blue lines). The photons are drawn from the SOURCE event class in 3.9 years of data, with a match radius of 0.15/0.3° for 
FRONT/BACK converting events. The expected ene rgies (red vertical lines) and histogram binning (black) arc identical to our previous paper 
on the G alactic center |Su fc Finkbcincr 2012b). The upper right panel is the same spectrum after convolving with the LAT line spread 
function jRaiaraman et al.ll2012l f and multiplying by E. The lower two panels are the same as the upper two panels, but for associated 
point sources. The associated source spectrum shows a bump between the blue lines due to the strong selection bias, but shows no evidence 
for a doublet at the expected energies. 



the LAT data processing as a source of the Galactic cen- 
ter lines. Even a lower significance detection would be 
interesting, because there is no trials factor for choice of 
line energies. 

Many of the Fermi-LAT point sources are associated 
with counterparts at other wavelengths, including blazars 
(BL Lacs, Flat Spectrum Radio Sources (FSRS), etc.), 
other AGNs (Seyferts, Radio Galaxies, e tc.), pulsars 
and binaries, and other Galactic sources ( Nolan et al.l 
Although su bstantially improve d over the First 



Fermi-LAT catalog ( Abdo et al. 2010ah . there are still 
344 sources in the 2FGL (—15% of the total) without 
obvious counterparts at Galactic latitude \b\ > 5°. 

Various statistical methods and theoretical scenar- 
ios have been suggested to classify and explain these 



unassociated sources (jAckermann et al.l l2012al) . includ- 
ing the existence of new types of source classes, e.g. 
dark matter subhalos. Numerical simulations suggest 
that within the Milky Way halo, dark matter subha- 
los form at all mass scales down to the simulation res- 
olution. Less massive halos might show themselves 
as gamma-ray so urces without significant emission at 



other wavelengths ( Belikov et al.ll201lUAckermann et al 
2012bHBucklev fc Hooperll2010l) . If such a signal were de- 



tected, it would be the first non-gravitational signature 
of dark matter. 

In this work, we use 3.9 years of LAT data to study the 
gamma-ray spectrum of the unassociated point sources in 
the 2FGL catalog to search for gamma-ray line emission. 
We find that the energy spectrum shows two lines at 111 
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Fig. 2. — Same as Figure[T] but using the ULTRACLEAN event class. 

GeV and 129 GeV. In Section [2] we describe our LAT 
data selection and analysis procedure. In Section [3] we 
examine the spectral line emission by various statistical 
tests and we summarize our main findings in Section |4] 

2. FERMI DATA SELECTION 

In this section, we briefly describe our data selec- 
tion and analysis procedure. We refer to our 



previ- 



ous papers for more deta iled information (|Su et al. 
Su fc Finkbeinenl2012alb[ ). The Fermi LAT is a pair- 



conversion telescope, in which incoming photons convert 
to e + e~ pairs, which are then tracked through the de- 
tector. The arrival direction and energy of each event 
are reconstructed, and the time of arrival recorded. The 
LAT is designed to survey the gamma-ray sky in the en- 
ergy range from about 20 MeV to several hundreds of 
GeV. The point spread function (PSF) is about 0.8° for 
68% containment at 1 GeV and decreases with energy 
as r 68 ~ E~ - 8 , asymptoting to ~ 0.2° at high energy. 
It is convenient to distinguish between front-converting 



and back-converting events that convert in the front and 
back regions of the tracker, respectively. The 68% con- 
tainment radius at high energy is r 68 ~ 0.15° for front- 
converting and res ~ 0.30° for back-converting events, 
with some dependence on the incidence angle on the de- 
tector. 

The LAT energy resolution (i.e. the half-width of the 
68% containment re gion) is of order 10% over mos t of 
the energy range ( see iFermi-L AT Collaborationll20 1 2i for 
details). Around 100 GeV, the resolution is closer to 7% 
for high incidence-angle events, and twice that for normal 
incidence. 

We use the latest publicly available data and instru- 
ment response functions, known as Pass 7 (P7_V6) 5 . We 
perform our analysis on both SOURCE and ULTRACLEAN 
events, and present figures based on each for compar- 
ison. The former has larger effective area and higher 

5 Details at http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
documentation/Pass7_usage .html 
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Fig. 3. — Upper panel: same as the upper left panel of Figure[l] 
(with an extra power of E, but for all unassociated sources. The 
green line is a power law dN/dE ~ _E -2,6 fit to events at E > 135 
GeV. Lower panel: same but for ULTRACLEAN class. The suspicious 
bump at logi?=2.2 in the upper panel disappears in the lower panel. 
However, the background estimate in the lower panel appears to be 
unreasonably low, so the green line is multiplied by a factor of 2.5 
to match the lower energy emission. Because of this disagreement, 
we do not use the ULTRACLEAN data for our final result. 



background. At the energies of interest CLEAN and 
ULTRACLEAN events give identical results. 

Photons coming from the bright limb at Earth's 
horizon, dominantly produced by grazing-incidence CR 
showers in the atmosphere, are a potential source of con- 
tamination. We remove this background by selecting 
events with zenith angle less than 105°. We also exclude 
some time intervals when data quality is poor, primarily 
while Fermi passes through the South Atlantic Anomaly. 



High resolution simulations suggest the existence of hi- 
erarchical dark matter halo structures on all r esolve d 

2008). 



mass scales 



e.g. 



Kuhlen ct al. 2009; Picri et al 



The high end of the mass function is visible in the 
Milky Way as dwarf galaxies, including the Magellanic 
clouds. Less massive subhalos could be too small to con- 
tain enough baryonic matter to be visible at other wave- 
lengths, but shine only via annihilation of dark matter 
particles in gamma rays. Such sources may appear as 
unassociated Fermi-LAT sources in the 2FGL catalog, 
which forms the basis for our line search. 

3.1. Source selection 

The Fermi-LAT 2FGL catalog consists of 1873 sources 
(100 MeV-100 GeV), of which 1290 are firmly identi- 
fied/associated an d 575 (31%) are unassociated sources 
( Nolan et af]l2012h . Cutting to Galactic latitude \b\ > 5° 
to avoid contamination from the Galactic disk results in 
344 unassociated sources. Since dark matter emission is 
expected to be non-variable in time, we also remove the 
25 among the 344 unassociated sources which have been 
flagged with variability_index > 41.6. Then we select 
photon events with energy in the range 100-140 GeV and 
zenith angle less than 105°. We include both FRONT and 
BACK converting events. Among the remaining 319 unas- 
sociated sources we select ones with at least one 100-140 
GeV photon within 0.15°/0.3° for FRONT/BACK convert- 
ing events, which results in 16 unassociated sources for 
SOURCE event class. The detailed information of these 16 
sour ces can be found in Table Q] 

In iNolan et al.l (|2012n . it has been noted that 51% of 
the unassociated sources have been flagged due to various 
issues (compared to 14% of the associated sources). We 
have tried cutting on each flag bit and find no significant 
impact on our results. In order to avoid introducing an 
additional trials factor, we decided not to cut on any 
flags. 



3.2. Composite energy spectrum 

We show the energy spectrum of the regions within 
0.15°/0.3° for FRONT/BACK converting events in two ways. 
In the left panels of Figure [1] and [2j we show the photon 
counts binned i n log E. The binning us ed is the same 
as in Fig. 16 of ISu fc Finkbeinerl (|2012bf ). This binning 
was chosen to optimize the signal on the 111 GeV and 
129 GeV lines in the Galactic center, and has not been 
modified for this study. In the right panels of Figure 
[U we plot E 2 dN/dE for the unbinned event di stribution 
convolved with the line-spread function (LSF) ( Edmondsl 
l201lL ISu fc Finkbeinerl l2012bT ). Each of these represen- 
tations has pros and cons. The histograms make it easy 
to see how many photons contribute to each bin, and al- 
low simple computations of Poisson likelihoods. The un- 
binned spectra give a sense of the spectrum without any 
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Fig. 4. — The inclination angle distribution of the 100-140 GeV photons in left panels of Figure [TJ and Figure [2] i.e. for SOURCE and 
ULTRACLEAN events, respectively. 
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Fig. 5. — The energy spectrum of unassociated 2FGL sources 
which have a 100-140 GeV photon within 0.15°/0.3° radius for 
FRONT/BACK LAT events. The spectrum is obtained from the 2FGL 
catalog. Sources thought to be potentially confused with Galactic 
diffuse emission are shown with dotted lines. The dashed black 
line shows the average of all the spectra. Each band shows integral 
photon flux from [0.1-0.3, 0.3-1, 1-3, 3-10, 10-100] GeV respectively 
from the likelihood analysis in that band with fixed photon power- 
law index. A 2a upper limit is shown if the source is not significant 
in a band. 

arbitrary binning choices, and convolution by the LSF 
allows maximum sensitivity to faint signals. However, 
the spectrum is "twice convolved" (once by the instru- 
ment and once by the processing) making the lines blend 
together to an undesirable degree. In the following we 
use the binned histograms for analysis, and provide the 



Fig. 6. — The angular distribution of 100-140 GeV photons 
matching with unassociated 2FGL sources within 0.15°/0.3° radius 
for FRONT/BACK LAT events (shown with solid/blue and dash/red 
curve, respectively) . The distribution is normalized by the annular 
area given a radius. The FRONT events shows more concentrated 
distribution than the BACK events, consistent with the expectation 
based on the point spread function. Both FRONT and BACK events 
suggest a central concentrated distribution. 

smoothed spectra only for reference. 

We plot the energy spectrum for photons near the 
16 unassociated sources with a photon at 100-140 GeV 
for SOURCE (Figured]) and ULTRACLEAN (Figured]) event 
classes. While this selection obviously suppresses the 
spectrum outside of the 100-140 GeV range, there is no 
way it can rearrange photons in the 100-140 GeV win- 
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Fig. 7. — Probability of obtaining the observed counts, in the 
energy bins centered on 111 and 129 GeV, in the Galactic center 
and subhalos as a function of the line fraction / = ^129/(^111 + 
-Fl29)- We find that the best fit ratio of the 129 GeV line to 111 
GeV line is 1.5, and the 2a range of the line ratio is [0.84, 4.5]. See 
Section 13.61 for details. 

dow. Intriguingly, we find two gamma-ray emission lines 
at 111 GeV and 129 GeV. One interpretation is that 
unassociated sources emit a gamma-ray doublet. An- 
other is that some flaw in the LAT data preferentially 
maps events to these energies. In order to test this hy- 
pothesis, we apply the same selection procedure to asso- 
ciated point sources, also shown in Figures Q] and [2j and 
find no line features at 111 GeV or 129 GeV. 

For each of the 16 unassociated sources, we show the 
integrated photon flux in energy bands (0.1-0.3, 0.3-1, 
1-3, 3-10, 10-100) GeV respectively in FigureO obtained 
from the 2FGL catalog. Among the 16 sources, three of 
them are marked with potentially confused with Galactic 
diffuse emission. 

3.3. Statistical significance 

Even with low statistics (7 counts at 111 GeV and 6 
counts at 129 GeV) it is possible to obtain a significant 
result if the backgrounds are low enough. Because WIMP 
annihilations can produce lower energy photons (final- 
state radiation, Z/W continuum, inverse-Compton, etc.) 
it may be incorrect to use lower energy emission to assess 
the background. However, at high energy there are very 
few photons in these sources, and there would be none 
from a 129 GeV WIMP. As a compromise, we assume 
the background is a power law, fit its amplitude to high 
energy (135 < E < 270), but choose the power-law index 
so that lower energy emission is modeled approximately 
correctly (Figure [3]). 

We assess the Poisson probability of observing 13 (or 
more) SOURCE counts in the two spectral bins with the 



background estimate in the upper panel of Figure[3] This 
has a probability of p — 0.00069 corresponding to 3.2tr. 
Removing sources to be potentially confused with Galac- 
tic diffuse emission (marked out in 2FGL) only mildly 
affect our results (3.3cr). The ULTRACLEAN events would 
give a much higher significance (> 4er) if we could believe 
the background estimate, but it looks implausibly low. 

3.4. Spatial distribution 

The subhalo candidates identified in this work are 
mostly distributed at |6| < 20°, at all longitudes. It 
is not clear whether this could be a selection effect, a 
fluke, or a hint about the true distribution of dark mat- 
ter subhalos. On one hand, dark matter subhalos pref- 
erentially dragged into the Galactic disk may lead to 
disk-l ike configurations, e.g. the propose d "dark disk" 



(e.g. lBruch et al.ll2009l; iPurcell et al.ll2009D . On the other 
hand, the distribution is not concentrated in longitude, 
so they may be nearby subhalos with lower mass, close 
enough to appear brighter than more massive subhalos, 
e.g. those hosting dwarf galaxies. 

3.5. Radial profile 

In Figure O we show the stacked angular distribution 
of 100-140 GeV photons, which are selected by match- 
ing with unassociated 2FGL sources within 0.15°/0.3° 
radius for FRONT/BACK LAT events, with respect to the 
source center provided by 2FGL. The distribution is nor- 
malized by the annular area at each radius. The FRONT 
events show a more concentrated distribution than the 
BACK events, consistent with the point spread function. 
Both FRONT and BACK events suggest a centrally concen- 
trated distribution. 

3.6. Line ratio 
Our previous work (|Su fc Finkbeinerll2012bh found 4 



(14) photons above background at 111 (129) GeV. This 
led us to expect the 129 GeV line might be stronger, 
but this work finds the 111 GeV to have slightly more 
counts: 6 (5) at 111 (129) GeV above background. Are 
these results compatible? 

In order to determine a confidence interval for the line 
ratio, we consider a total of N photons for the doublet, 
with k of them in the 129 GeV bin, and the rest in the 
111 GeV bin. The binomial probability of observing k of 
N counts in this bin is 

where / = -£129/(^111 + -F129) is the true fraction of dou- 
blet photons at 129 GeV. Figure [7] shows this probability 
(i.e., the probability of observing k counts given N and 
/) as a function of / for the GC, subhalos, and the prod- 
uct of the two. 
To obtain a confidence interval, we find /i ow such that 

N 

P(k > x, n, / low ) = Pb(k, n, /) = 0.025 (2) 
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Fig. 8. — The spatial distribution of unassociated 2FGL sources which have a 100-140 GeV photon with in 0.15° /0. 3° radiu s for FRONT/BACK 
LAT events (red stars). Fo r comparison, we also show t he spatial distribution of 16 dwarf galaxies (e.g. lAbdo et al.ll2010bl . blue triangles), 
106 nearby galaxy clusters (Rciprich & Bohrinecr 2002, dark green circles), and associated sources with the same matching criteria (purple 
squares). We find no spatial overlaps between these sources. Unassociated sources thought to be potentially confused with Galactic diffuse 
emission (red circles) and unassociated sources with less than 5<r detection (red squares) are noted. The lower panel shows the same plot as 
the upper panel but selected unassociated point sources by matching with 150-500 GeV events. In both cases, the sources with high-energy 
photons are mostly near the disk. This implies that the Galactic latitude distribution may result from a selection effect. 



with a complementary expression for /high- The 95% 
confidence interval (corresponding to "2cr" confidence) is 
then 0.167 < / < 0.765 for the subhalos and 0.524 < 
/ < 0.935 for the Galactic center. A significant range of 
/ is allowed in both cases, so we can combine the counts 
from both and obtain 0.457 < / < 0.820 for the joint 
fit. This yields 95% confidence bounds on the line ratio 
0.84 < F129/.F111 < 4.5. The data are consistent (at 2a) 
with the lines being equally strong, but also with the 



129 GeV line being 4.5 times as strong. Clearly more 
data will be required to measure the line ratio with high 
confidence. 

4. DISCUSSION AND CONCLUSION 

In this paper, we have reported evidence for line emis- 
sion at 111 GeV and 129 GeV from unassociated Fermi- 
LAT point sources. The lines have a significance of 
p = 6.9 x 10~ 4 or 3.2cr for a simple power-law background 



model. These results provide independent support for 
our previous claims of a double gam ma-ray line in the 
Galact ic center at the same energies (|Su &: Finkbeiner 
2012bl) . The double line emission is compatible with the 
scenario of a 129 GeV WIMP annihilating to 77 and 
7Z, producing the two lines. As a test of systematics, 
we apply the same selection and analysis procedure to 
associated Fermi-LAT point sources, and find no evi- 
dence for lines at these energies. It is difficult to imagine 
instrumental systematics that could produce this double 
line emission only at the Galactic center region and the 
locations of unassociated point sources without affecting 
other regions of the sky. We find this evidence to be 
persuasive, but it cannot be considered conclusive until 
more data become available. 

Further observations are essential, not only to firmly 
establish the existence of the lines, but to measure the 
line ratio. The ratio of jZ and 77 line strength depends 
only on the particle physics model, and is independent 
of the astrophysical uncertainties such as the dark mat- 
ter distribution in the halo. The line pair is also com- 
patible with a 141 GeV WIMP annihilating through 7Z 
and 7/1 f or mn ~ 125 GeV, a s in the "Higgs in Space" 



scenario (Jackson et al 



2010). 



However, 



we have not 
141 GeV. In 



found any significant gamma-ray line at 
any case, additional data will be critical for measuring 
the line ratio, which is currently only poorly determined 
(Figure El). 

A possible change to the Fermi scan strategy could ac- 
cumulate S/N on the double spectral lines in the Galac- 
tic center up to ^4 times as fast as the current survey 
strategy, and it is crucial for studying the double line 
emission. We believe this could be done with only mod- 
est impact on other Fermi science objectives. With the 
huge effective area and low energy threshold of H.E.S.S 
II, it may be possible to confirm a spectral bump in the 



Galactic center fairly soon. However, the energy resolu- 
tion of H.E.S.S. II is inferior to LAT at 129 GeV, and it 
may be difficult to resolve the doublet. 

The 2FGL catalog is based on 24 months of LAT data, 
and an updated catalog based on 48 months would pro- 
vide improved source parameters and associations, de- 
creasing the background noise for the subhalo analysis 
presented here. Furthermore, multi-wavelength follow- 
up observations would be helpful to identify the nature 
of unassociated 2FGL sources and refine the list of asso- 
ciations. 

By stacking the unassociated 2FGL point sources to- 
gether, it might be possible to reveal the spatial profile 
of the gamma-ray distribution at 111 GeV and 129 GeV. 
One may in principle improve the positional informa- 
tion of these unassociated sources by using high energy 
photons and better quantifying the gamma-ray spatial 
profile. However, with ~ 1 photon per source at high 
energy, the details of the algorithm used for centroiding 
the sources are critical, and consideration of the spatial 
profile is beyond the scope of this work. We simply use 
the position provided by the Fermi 2FGL catalog. 
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TABLE 1 



Source name 


T? A 





£ 


v. 
D 


r lelgS 


(7 


Variability 


Spec index 


txadius 


Spec type 


2FGL 


J0341.8+3148c 


55.5 


31.8 


160.3 


-18.4 


32 


7.6 


24.7 


2.31 


± 





09 


0.0706 


PowerLaw 


2FGL 


J0526.6+2248 


81.7 


22.8 


182.9 


-6.9 


1 


7.1 


27.5 


2.88 


± 





16 


0.0641 


PowcrLaw 


2FGL 


J0555.9-4348 


89.0 


-43.8 


250.4 


-28.0 





5.0 


29.1 


2.11 


± 





17 


0.0994 


PowerLaw 


2FGL 


J0600.9+3839 


90.2 


38.7 


173.2 


7.6 





5.1 


14.1 


2.04 


± 





21 


0.0490 


PowerLaw 


2FGL 


J1240.6-7151 


190.2 


-71.9 


302.1 


-9.0 





8.2 


18.4 


1.82 


± 





16 


0.0458 


PowerLaw 


2FGL 


J1324.4-5411 


201.1 


-54.2 


307.8 


8.4 


24 


4.9 


23.3 


2.39 


± 





14 


0.1126 


PowerLaw 


2FGL 


J1335.3-4058 


203.8 


-41.0 


311.8 


21.1 





4.6 


11.2 


2.13 


± 





18 


0.0800 


PowcrLaw 


2FGL 


J1601. 1-4220 


240.3 


-42.3 


336.3 


7.9 





7.3 


20.8 


2.46 


± 





10 


0.1037 


PowerLaw 


2FGL 


J1639.7-5504 


249.9 


-55.1 


331.8 


-5.6 


9 


5.9 


21.1 


2.79 


± 





14 


0.0568 


PowerLaw 


2FGL 


J1716.6-0526c 


259.2 


-5.4 


16.6 


18.2 


2080 


6.8 


25.7 


2.43 


± 





26 


0.1052 


LogParabola 


2FGL 


J1721.5-0718c 


260.4 


-7.3 


15.6 


16.2 


41 


6.0 


20.9 


2.68 


± 





36 


0.0795 


LogParabola 


2FGL 


J1730.8+5427 


262.7 


54.5 


82.2 


33.3 





4.6 


16.9 


2.69 


± 





18 


0.1258 


PowerLaw 


2FGL 


J1844.3+1548 


281.1 


15.8 


46.3 


8.7 


4 


12.5 


29.8 


2.43 


± 





08 


0.0403 


PowerLaw 


2FGL 


J2004.6+7004 


301.2 


70.1 


102.9 


19.5 





9.5 


36.8 


1.97 


± 





11 


0.0368 


PowerLaw 


2FGL 


J2115.4+1213 


318.9 


12.2 


62.6 


-24.5 





5.1 


25.3 


2.38 


± 





19 


0.0800 


PowerLaw 


2FGL 


J2351.6-7558 


357.9 


-76.0 


307.7 


-40.6 





4.1 


20.8 


1.92 


± 





19 


0.0702 


PowcrLaw 



Note. — The table provides detailed information about the unassociated 2FGL point sources we have identified with at 
least one 100-140 GeV photon within 0.15/0.3° for FRONT/BACK events. The first column is the 2FGL catalog name in the 
format 2FGL JHHMM . m+DDMM, where V indicates that the source is considered to be potentially confused with Galactic diffuse 
emission. The second/third column are Right Ascension (J2000) and Declination (J2000). The forth and fifth column are 
Galactic Longitude and Galactic Latitude. The sixth column is the flag parameter which indicate possible issues noted in 
detection or characterization of the source. Sources having no flags raised with value are those without potential problems. 
The seventh column is the variability index, defined as the test statistic for the hypothesis that monthly averages of the 
source flux vary relative to the null hypothesis of constant flux. The TS is distributed as x 2 with 23 degrees of freedom, so 
a value greater than 41.64 indicates a > 99% chance of being a variable source, and we have removed these sources from our 
analysis. The eighth column shows the best fit for the photon number power-law index (for logarithmic parabola spectra 
it is index at the Pivot Energy) derived from the likelihood analysis for 100 MeV-100 GeV. The ninth column shows the 
average of semimajor/semiminor axis of the error ellipse at 68% confidence. Source detection significance in Gaussian cr 
units is shown in the tenth column, which is derived from the likelihood Test Statistic for 100 MeV-100 GeV analysis. The 
eleventh column shows the best fit form of the spectral type. We note that only three sources have logarithmic parabolic 
spectral shape (two of them are marked with potenti ally confused with Galactic diffuse emission). Detailed explanation of 
parameters listed in this table can be found in lNolan et al. | J2012T) . 

APPENDIX 

We show the same energy spectrum as in Figure Q] but with high incidence angle events only with 9 > 40° in Figure 
E3 In FigureUm we show that the 111 GeV and 129 GeV lines do not appear in the SOURCE minus ULTRACLEAN events, 
i.e. cosmic ray contamination is not a plausible explanation for the line emission. In Figure [Til and Figure [T^] we show 
the unassociated/associated point sources by matching with LAT events of different energy range. We also compare 
with distributions of all the unassociated/associated sources. Selection effect is plausible explanation for the spatial 
distribution of selected unassociated sources shown in Figure [H 
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Spectrum of unassociated 2FGL sources 



: Eselect: 100-140 GeV 



c 3 



1 - 




1.2 



1.4 



1.6 1.8 2.0 

log(E/GeV) 



2.2 



2.4 



Spectrum of unassociated 2FGL sources 



: Eselect: 100-140 GeV 



3. 2 



1.2 



1.4 
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log(E/GeV) 



2.2 2.4 



Fig. 9. — Same as the upper left panel of Figure [T] but using only events with high incidence angle 9 > 40° which has better energy 
resolution to reveal the energy spectrum of unassociated 2FGL catalog. We found two gamma-ray line emission on 111 GeV and 129 GeV. 




Fig. 10. — Same as the upper left panel of Figure[T] but using only events belongs to SOURCE class but not to ULTRACLEAN event. This 
selected set of events are dominated by cosmic ray events. We found no significant gamma-ray line emission on 111 GeV and 129 GeV 
(0.47it with background estimation using 100-500 GeV data). 



11 



Eiielect: 150-500 GeV 



0.4 0.6 
sin(b) Galactic latitude 




E 

z 2 



Eselect: 60-100 GeV 
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Fig. 11. — The distribution of unassociated 2FGL sources selected by matching with LAT photons of different energy range. The panels 
from the upper left to lower right are for 150-500 GeV, 100-140 GeV, 60-100 GeV, 30-60 GeV, and 10-30 GeV, and for comparison all 
unassociated sources with |b| > 5° and variability index < 41.64. 
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Fig. 12.— 



The same as Figure [TT] but for the associated 2FGL sources. 



